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1: 

■ Introduction 

> : 

In this paper, we give an algorithm to compute the following cohomology 
t-h ■ groups on U — C n \ V(f) for any non-zero polynomial / G Q[x±, . . . , x n ]; 

1. H k {U, Cjj), C[/ is the constant sheaf on U with stalk C. 

On i 

We also give partial results on computation of cohomology groups on U for a 
locally constant sheaf of general rank and on computation of H k (C n \ Z, C) 
>• ! where Z is a general algebraic set. 

Our algorithm is based on computations of Grobner bases in the ring of 
differential operators with polynomial coefficients, algorithms for functors in 
the theory of P-modules ( [p4}| and [[^|), and Grothendieck-Deligne compar- 



2. H k (U, V), V is a locally constant sheaf of rank 1 on U . 



ison theorem |12[, PJ, which relates sheaf cohomology groups and de Rham 
cohomology groups. 

One advantage of the use of the ring of differential operators in algebraic 
geometry is that, for example, Q[x,l/x], which is the localized module of 
Q[x] along x, is not finitely generated as a Q[x]-module, but it can be re- 
garded as a finitely generated Q(x, cQ-module where d x = d/dx. In fact, we 

(\ k — 1 
~. Computation of the localization 
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of a given D-module and computation of the integration functor are the most 
important part of our algorithm. 

As an introduction to this paper, it will be the best to mention about 
how we started this project. 

Let us consider the differential equation for the function / = (x — u) a (x — 
v) b where u and v are rational numbers (u < v) and a and b are rational 
numbers such that a+b (jL Z. The function / satisfies the differential equation 

pf — 0, p = (x — u){x — v)d x — a(x — v ) — b(x — u). 

Let p be the formal Fourier transform of this operator 

P = {-d x - u) (-d x - v)x - a(-d x - v) - b(-d x - u) 

= xd 2 x + (ux + v x + 2 + a + b)d x + uvx + u + v + av + bu 

and A be the ring of differential operators Q(x, d x ). We want to evaluate the 
dimension of the Q- vector space 

A/(Ap + xA) ~ (A/Ap)/x(A/Ap) ~ (A/Ap)/lmx-, 

which is called the (0-th) restriction of A/Ap along x = 0. Now, we can 
apply the algorithm for the D-module theoretic restriction in Section 5] 
to evaluate the dimension. Here, we need what is called a 6-function for the 
evaluation, and it is nothing but the indicial (characteristic) polynomial at 
x = of the ordinary differential operator p that appears in the classical 
method of Frobenius; the polynomial is s(s — a — b). Applying Proposition 
5.2 in with this 6-function, we conclude that the dimension is equal to 
one, which is the number of the bounded segments of R \ {u, v}. 

Next, we tried to evaluate the dimension of A 2 / (A 2 p + A 2 q + xA 2 + yA 2 ) 
where A 2 is the ring of differential operators generated by x,y,d x and d y , 
and p and q are differential operators which annihilate the function / = 
x a y b {l — x — y) c ; p = x(l — x — y){d x — a/x + c/(l — x — y)) and q = 
y(l — x — y)(d y — b/y + c/(l — x — y)). We evaluate the dimension, this 
time with a computer program ||33|| , by iterating to apply the algorithm for 
computing the 0-th restriction in pi] firstly for x and secondly for y. The 
result is again one, which is equal to the number of the bounded cells of the 
hyperplane arrangement R 2 \ {(x,y) \ xy(l — x — y) = 0}. It is well known 
in the theory of hypergeometric functions that the number of bounded cells 
is equal to the dimension of the middle dimensional twisted cohomology 
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group associated with /, which is equal to the rank of the corresponding 
hypergeometric system. (Strictly speaking, it turns out that p and q do not 
generate the annihilating ideal {£ £ A2 \ if = 0} for / (see Example 
however, the ideal generated by p, q happens to be 'close enough' to the 
annihilating ideal.) 

Inspired by the observation above, we started the project to obtain an 
algorithm for computing the cohomology groups of the complement of an 
affine variety by elaborating the method sketched above. 

1 Computation of cohomology groups of the 
complement of an affine hypersurface 

For any non-zero polynomial f(x) £ Q[si, . . . ,x n ], we prove the following 
theorem. 

Theorem 1.1 Put X = C n , Y = V(f) := {x £ X | f(x) = 0}, and 
U = X \ Y . Then the cohomology group H k (U,Cu) is computable for any 
integer k, where Cu denotes the constant sheaf on U with stalk C. 

Note that H k (U,Cu) is the k-th cohomology group (with coefficients in 
C) of the 2n-dimensional real C°°-manifold U c i underlying U. Also note that 
H k (U, Cu) = for k > n since U is affine and that H°(U, C v ) = C since U 
is connected. 

In the theorem above, we may replace Q by any computable field. Here, 
we mean by a computable field a subfield K of C such that each element of 
K can be expressed by a finite set of data so that we can decide whether 
two such expressions correspond to the same element, and that the addition, 
subtraction, multiplication, and division in K are computable by the Turing 
machine. For example, any algebraic extension field of Q of finite rank is 
a computable field by virtue of Grobner bases and factorization algorithms 
over algebraic number fields. 

In this section, we illustrate an algorithm to compute the cohomology 
groups. Correctness will be proved as a special case of the corresponding 
theorem for cohomology groups with the coefficients in a locally constant 
sheaf of rank one. In order to compute the cohomology groups, we translate 
the problem to that of computations of functors, especially to that of the 
de Rham functor, of ^-modules, which are studied in a series of papers 
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23 . Here, A n is the ring of differential operators with polynomial 

coefficients and is called the Weyl algebra. The computations of functors 
are based on the Buchberger algorithm to compute Grobner bases in the 
Weyl algebra. We shall quickly review the definition of Weyl algebra and the 
Grobner basis. See PHI , ||, ITfl , f2lfl for details, |30|] for an introduction, 
and |33] for implementations. 
The Weyl algebra 



A n — Q(xi, . . . , x n , di, . . . , d n ) 

is the ring of non-commutative polynomials generated by 2n elements Xi,di, 
(i — 1, . . . , n) satisfying the relations 



OC 2 30 j 



1 (i = j) 




XjOi 



The theory (and practice) of Grobner bases works perfectly well for left 
ideals in the Weyl algebra A n . We quickly review the relevant basics. Every 
element p in A n can be written uniquely as a Q-linear combination of nor- 
mally ordered monomials x a d b . This representation of p is called the normally 
ordered representation. For example, the monomial d\Xidi is not normally 
ordered. Its normally ordered representation is Xidf + d\. 

Consider the commutative polynomial ring in In variables 

gr(A„) = Q[xi,...,a; n ,fi,...,f n ] 

and the Q-linear map gr : A n — > gr (A n ), x a d h t— * x a C, b . Let < be any term 
order on gr (A n ). This gives a total order among normally ordered monomials 
in A n via x A d B > x a d b x A ^ B > x a £, b . For any element p G A n let in < (p) 
denote the highest monomial x A d B in the normally ordered representation 
of p. If J is a left ideal in A n then its initial ideal is the ideal gr (m<(J)) in 
gr (An) generated by all monomials gr (in < (p)) for p E I. Clearly, gr (m<(J)) 
is generated by finitely many monomials x a ^ b . A finite subset G of / is called 
a Grobner basis of / with respect to the term order < if {gr (in < (q)) | q E 
G} generates gr(m<(/)). Noting that in < (p) < in < (q) implies in < (hp) < 
in < (hq) for all h E A n , one proves that the reduced Grobner basis of I is 
unique and finite, and can be computed using Buchberger's algorithm. Any 
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left (or right) ideal in A n is finitely generated and we denote by (pi, . . . ,p m ) 
the left ideal in A n generated by pi, . . . ,p m G A n . 

Most constructions in the commutative algebra can be reduced to com- 
putations of Grobner bases. This is also the case with some constructions 
for modules over the Weyl algebra. For example, the construction of a free 
resolution of a left coherent A n -module is a straightforward generalization of 
algorithms of constructing free resolutions of modules over the ring of polyno- 
mials. As to algorithms to construct a free resolution by the Schreyer order, 
see |], p. 167 Theorem 3.7.13], || Theorem 15.10] and f29|. We note that an 
algorithm to compute a sheaf cohomology on the n-dimensional projective 
space is given based on computation of syzygies in the ring of polynomials 
([PH). Computation of an elimination ideal in the Weyl algebra is also a 
straightforward generalization of computation of an elimination ideal in the 
ring of polynomials, (see, e.g., [TIL p. 69 Theorem 2.3.4], J7L p. 114 Theorem 2]). 
These two constructions will be used in our algorithm to obtain cohomology 



groups (see Algorithm [L2| Step 3, Procedure Step 2 and Procedure |272 
Step 2). 

However, the non-commutativity causes some difficulty in constructing 
various objects in the category of modules over the Weyl algebra. For ex- 
ample, to compute the tensor product of right and left A n -modules in the 
derived category, special care must be taken. This problem has been an open 
problem since p2| . As a special (but important) case of the tensor product 
computation as above, we give an algorithm for the D-module theoretic re- 
striction of an A n -modu\e by using the F-filtration and the 6-function (or the 
indicial polynomial). As to details, see []24], Proposition 5.2, Theorem 5.7, 



Algorithm 5.10] and |[25|| . U.Walther ]34|] solved a related problem of com- 
puting algebraic local cohomology groups based on ^/-filtration, 6-function 
and the Cech complex. As we will see in Section [7|, his algorithm gives an 
algorithm for Theorem |1.1| different from ours explained below. 

We have explained a general background on an algorithmic treatment of 
modules over the Weyl algebra. Now let us explain our algorithm to compute 
cohomology groups by a top-down expansion. 

Algorithm 1.2 (Computation of the cohomology groups H k (U, Cu)) 
Input : a polynomial / G Q[ 

Output : H k (U, Cu) for < k < n where U = C n \ V(f). 
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1. Find a left ideal I such that 

Q[x, j)^A n /I 

as a left y4 n -module. 

2. Let J be the formal Fourier transform of /; 

J = h 

3. Compute a free resolution of length n + 1 

A P-(n+l) -L-^+V A p_ n A P-(n-l) -L- 1 A p , I r p 

^ln ' -ri n ' Sin ' Sl n ' -f±n/ " ' u ) 

(po = 1) of A n /J by using Schreyer's theorem [|, Theorem 15.10] with 
an order which refines the partial order defined by the weight vector 

di ■■■ d n xi ■■■ x n 
( 1 ■•■ 1 -1 ••■ -1 ). 

4. Compute the cohomology groups of the complex of Q-vector spaces 

(A n /( Xl A n + ■■■ + x n A n ) ® An AP n ~\ 

Then, the (k — n)-th cohomology group Ker (1 £g>L fc_n )/Im (1 (g)L fe ~ n ~ 1 ) 
of the complex above tensored with C gives H k (U, Cy). 



The step 1 will be explained in Procedure 1.3 in detail and the steps 2, 3 



and 4 will be explained in Procedure |1.4] in detail. 



We note that the steps 2, 3 and 4 are nothing but the computation of 
H k - n (A n /(d 1 A n + ■■■ + d n A n ) ® L An A n /I), 
which is denoted by 

rk—n fk—n 

/ A n /I= Q[x,l/f] 



in the theory of ©-modules. We shall prove that this cohomology group ten- 
sored with C is equal to H h (U, Cjj) by the Grothendieck-Deligne comparison 
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theorem in Section |^. Here, for a left y4„-module M and a right A n -module 
N, we denote by N ®\ n M the complex 

(N ® An M\ 1 <g> d*" 1 ; 2 = 1, 0, -1, -2, . . .) 

where 

• • • M i ■ ■ ■ M' 1 — >M — ► 0, (exact) 



M % is a free y4„-module, d % is a left A n -morphism and M 1 = 0, d° = 0. It 
is known that there exists a finite length free resolution for a given finitely 
generated left A n -module M (e.g., apply the method of || p. 336 Corollary 
15.11] to our case). 

Remark 1.1 For a left v4 n -module M = A n /I, the left A„_!-module (A n /d n A r 
M = A n /(I + d n A n ) is called the 0-th integral of M with respect to x n . Why 
is it called the integral? Let us explain an intuitive meaning of this termi- 
nology. 

Let / be a function of xi,...,x n . We suppose that the function / is 
rapidly decreasing with respect to the variable x n and put I = Ann f = {£ E 
A n \£f = 0}. Then the A n -module generated by the function / is isomorphic 
to the left y4 n -module A n /I. Put g(xi, . . . , x n -i) — f™oc fi x ii ■ ■ ■ > x n)dx n . 
Then, we have [(I + d n A n ) fl A n ^i]g = 0. In fact, since any element £ in 
(J + d„A n ) n v4 n _i can be written as £ = l\ + <9 n £2, ^i £ ^,-^2 G A n , we 
have = JZoi^i + d n £ 2 )fdx n = d n (£ 2 f)dx n = 0. Therefore, g can 
be regarded as a solution of the differential equations corresponding to the 
left v4 n _i-submodule A„_i/(I + d n A n ) n of A„/(d n A n ) ® An M. Note 
that An/ (d n A n ) ®A n M itself describes a system of differential equations for 
I-oo x if dx n with j > 0. 



Let us explain in detail the step 1 of Algorithm |1.2| . This algorithm is 
given in |2"3( 



Procedure 1.3 (Computing the differential equations for 1// r °; step 1 of 
Algorithm |1.2|) . 



Input: /. 

Output: a left ideal / of A n such that Q[x, 1/ f) ~ A n /I. 
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1. (Computation of the annihilating ideal of f s ) 
Compute 

(t - f(x), ^-d t + d l ,..., ^Ldt + d n ) n Q[td t ](x, d x ). 

Replacing td t by — s — 1, we obtain the left ideal Ann/ S in Q[s](x, d x ). 
(Call Procedure |4.1| with d — 1 to compute the intersection of the left 
ideal and the subring Q[td t ](x,d x ).) 

2. (Computation of the 6-function of /) 
Compute the generator b(s) of 

(Ajmf',f)nQ[s] 
by an elimination order x, d x > s. 

3. Let r be the minimum integral root of b(s) = 0. Put / = (Ann/ s ) s ^ rQ . 
Then, we have Q[x, j] ~ A n /I. 

The polynomial b(s) is called the (global) Bernstein- Sato polynomial or 
the b-f unction of /, which coincides with the (global) indicial polynomial of 
5(t — f), and plays an important role in our algorithms. This polynomial is 
the minimal degree polynomial satisfying the relation 

Lf s+1 = b(s)f, 3L G Q[s](x,d x ). 

The left module A n /I is a holonomic A ra -module (or called a module belong- 
ing to the Bernstein class). The holonomicity of A n /I and the existence of 
the 6-function were shown by I.N.Bernstein. See, e.g., [ffl] and ]|, p. 13, 5.5 
Theorem]. It is known that when / ^ const, b(s) always has s + 1 as a fac- 
tor. M.Kashiwara proved that all the roots of b(s) = are negative rational 



numbers for any f & C[x\, . . . , x n ] |T5 



Remark 1.2 The left y4 n -isomorphism Q[x, 1/f] A n /I is expressed as 
ip = ip 2 o tp± where the left A n -isomorphisms ipi and if2 



Q[ Xj i/f)j^A n r°^A n /i 

are defined as 

MD = i e A n /I, MD = f ro e KP 
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and 



f r -fcN _ L(r - k) ■ ■ ■ L(r - 1) 



Hence, for example, we have 

^(^) = roidf/dx^r- 1 = <pZ\r {df/d Xi )L{rv - l)r°/^o " !))■ 

Example 1.1 For / = x(l — x), we have 

Ann / s = (x(l - x)^ - s(l - 2a:)). 
The 6-function of / is s + 1 with 

((l-2x)d x + 4(l + s))f s+1 = (s + l)f, 

and hence we get 

Q[x, 1/f] ~ Q(x, <9 ffi )/(x(l - x)^ + (1 - 2z)). 

Example 1.2 Put / = x 3 — y 2 . We compute the left ideal I such that 
Q[x,y, 1/f] — Aij I. Here is a log of the output of kan/kO, which may be 
self-explanatory. The system kO is a translator that compiles Java like inputs 
to codes for kan/sml, which is a Postscript like language for computations 
in the ring of differential operators 



In(9)= a = annf s(x~3-y~2, [x,y] ) ; 
Computing the Groebner basis of 

[ v*t+x~3-y~2 , -v*u+l , -3*u*x~2*Dt+Dx , 2*u*y*Dt+Dy ] 
with the order u, v > other elements. 



In(10)=a : 

[3*x"2*Dy+2*y*Dx , -6* (-1-s) -2*x*Dx-3*y*Dy-6] 

In ( 1 1) =b=ReducedBase (Eliminatev (Groebner (Append [a , y~2-x~3] ) , 

[x,y,Dx,Dy])); 

In(12)= b: 

[ -216*s~3-648*s"2-642*s-210 ] 
In(13)=Factor(b[0] ) : 

[[-6,1], [6*s+5 , 1 ] , [6*s+7 , 1 ] , [s+1 , 1 ] ] 
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Since s = — 1 is the minimum integral root of the 6-function, we have 
Q[x, y, 1/f] ~ A 2 /(3x 2 d y + 2yd x , -2xd x - 3yd y - 6). 
Finally, let us explain our algorithm for computing 

H k - n (A n /(d 1 A n + ■■■ + 8 n A n ) ® L An AJI). 
This is a detailed explanation of steps 2, 3 and 4 of Algorithm p~^. We can 



compute the cohomology groups by applying [ 25| , Theorem 5.3] to the Fourier 
transformed ideal J of I. Correctness will be discussed in Sections [| and |. 
Put 

di ■■■ d n xi ■■■ x n 
w = ( 1 ••• 1 -1 ••• -1 ). 



and 
where 



F k = {/ G A n | ord w (f) < k} 



ord w (x a d b ) = —\a\ + 
{-Ft} is called the V^-filtration. 

Procedure 1.4 |^5[ (Computing the D-module theoretic integral of A n /I; 

steps 2, 3 and 4 in Algorithm |1.2j ) 

Input: a left ideal / of A n . (A/ 1 is holonomic.) 

Output: The — k-th cohomology groups of A n / (diA n + • — h d n A n ) A n /I 
for < k < n. 

1. Let J be the formal Fourier transform of /; 

^ ^\x i ^-8 i ,a i ^x i ,(i=l,...,n)' 

2. Let G be a Grobner basis of the left ideal J with the weight vector w. 
Find the generator b{9\ + ■ • • + 9 n ) of 

(in w (G))nQ[6 1 + --- + 9 n ], 6 i =x i d i . 

3. Let k\ be the maximum integral root of b(s) = 0. If there exists no 
integral root, then quit; the cohomology groups are all zero in that 
case. 
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4. Let < w be a refinement of the partial order by w. Construct a free 
resolution 

A P-(n+i) -L-^+V A p_ n -L^l 4 P_(„_i) -L- 1 4 p , / T n 

with p Q = by using the Schreyer orders associated with < w . 

5. (Computation of degree shifts) Put s\ = and 



i<?<p 



max (ord^L^ 1 ') + sf) (1 < % < p-(k+i)) 



successively. 

6. Compute the cohomology groups of the induced complex 

•■■ ^ F kl _ sl /(F^ + xA n )®---@F kl _ sl J(F^ + xA n ) 
^ F fcl /(F_ 1 + xA n )^0 

as a complex of Q- vector space where xA n = x\A n + • ■ • + x n A n . Then, 
the (k — n)-th cohomology group 

Keil k - n /\mL k ~ n - 1 

of this complex tensored by C gives H k (U, Cu)- 

In step 2, we denote by in w (J2( a ,b)ei c a bX a d b ) the w-leading form 

V" c a bX a d b , m = max (w, (a, b)) 

(w,(a,b))=m 

where (•, •) denotes the standard inner product in Z 2n . Put in w (G) = 
{in w (g) \g G G}. One needs to compute the intersection of the left ideal 
in w (J) = (in w (G)) and the subring. This can be done in a procedure similar 
to the one explained in Section [|. 

In step 4, we compute Grobner bases with the Schreyer orders over the 
order < w to construct a resolution. Note that < w is not a well-order, which 
causes a difficulty of computation of Grobner basis in our non-commutative 
situation. There are two ways to overcome this difficulty; one is to use the 
F-homogenization introduced in |2(| (see also p2| , Section 3]) and the other 
is the use of the homogenized Weyl algebra which has a homogenization 



11 



variable h so that the relation <9j£j = Xidi + h 2 holds. The homogenized Weyl 
algebra was introduced in kan/sml [33] since version 2 released in 1994. See 
on a theoretical study on this homogenization technique. 

In Step 6, we truncate the complex from above by using k\] we could also 
make truncation from below by using the minimum integral root of b(s) = 0, 
which would somewhat reduce the complexity of Step 6. See |25[] for details. 
Moreover, if we need to compute H k (U, Cu) only for k > £, then a resolusion 
of length n — £ + 1 suffices in Step 4. 

Example 1.3 We take / = x(l — x) in Q[x]. We denote by A the Weyl 
algebra Q(x, d x ). As we have seen, we have Q[x, 1/ /] ~ A/ (p), p = x(l — 
x)d x — (2x — 1). The formal Fourier transform of p is p = —xd 2 — xd x . By 
multiplying — x from the left, we have 

-xp = 6{6 - 1) + x6, 6 = xd x . 

Therefore the 6-function is equal to s(s — 1) and k\ = 1. The resolution of 
A/(p) is 

— > A - {x9 l^ dx) A — ► A/(p) — > 0. 

Since k\ = 1 and the degree shift by xd^ — xd x is equal to 1, the truncated 
complex is 

— > F /(F_! + xA) - {xd tf^ F l /{F_ l + xA) — > 0. 

Since 

F /(F_! + xA) = Q, + xA) = Q + CR 

and 1 • — xd x ) = in F 1 /(F_i + xA), we conclude that 

tf- 1 = Fo/(F_! + xA) = Q, F° = F /(F_! + xA) = Q 2 . 

Hence, the cohomology groups of U = C \ {0, 1} are 

H°{U,Cu) = C, H\U,Cu) = C 2 . 

The two generators of H 1 correspond to two loops that encircle the points 
x = and x = 1 respectively in view of the Poincare duality of homology 
groups and cohomology groups. 
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Example 1.4 (Cohomology groups of C 2 \ V(x 3 — y 2 )) 
This is the output of kan/kO. 

In(43)= bb=bf unctionFor Integral ( [3*x~2*Dy+2*y*Dx,-2*x*Dx-3*y*Dy-6] , [x,y]); 
In (44)= bb: 

[ -216*s~3+432*s~2-264*s+48 ] 
In (45)= Factor (bb): 

[[-24 , 1 ], [ 3*s-2 , 1 ], [ s-1 , 1 ], [ 3*s-l , 1 ]] 

In (46) =integralOf Module ( [3*x~2*Dy+2*y*Dx , -2*x*Dx-3*y*Dy-6] , [x,y] ,1,1,2) : 

Here, 1,1,2 specify the mimimum and the maximum integral roots, and the 
length of the resolution respectively. 

0-th cohomology: [ , [ ] ] 
-1-th cohomology: [ 1 , [ ] ] 
-2-th cohomology: [ 1 , [ ] ] 

The output means that 



H°(U,C U ) = C, H\U,C V ) = C, H 2 (U,C 



u 



0. 



y 2 , we 



Let us explain this example a little more precisely. For f = x 3 
have Q[x,y, 1/ f] ~ A 2 /I with 

/ = (2xd x + 3yd y + 6, 3x 2 d y + 2yd x ). 

Its Fourier transform is A 2 j J with 

J = (-2xd x - 3yd y + 1, 3yd 2 x - 2xd y ). 

The 6-function of A 2 / J is (s — I) (3s — 1) (3s — 2). Hence we put ki — 1. The 
free resolution of A 2 /J is given by 







A 2 -> A 2 /J -> 



with 



—2xd x — 3yd y + 1 
3yd 2 x - 2xd y 
-9y 2 d y d x - 3yd x - 4x 2 d y 
-27y 3 d 2 - 27y 2 d y + 3y + 8x 3 d y 





/ 3yd x 


2x 


-1 


\ 


L- 2 = 


2xd y 


-3yd y + 2 


-dx 







-9y 2 d y - 3y 





2x 


1 




V AX 2 8y 





-3yd y + 4 


d x j 


L- 3 = ( 


-3yd y + 2, 


-2x, d x , 


1 )• 
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The shift vectors are given by 

(s{, 4 si, s\) = (0,1,0,-1), 
(si, si, si, si) = (0,1,-1,0), 
si = 0. 

By computing the truncated complex, which is a complex of finite dimen- 
sional vector spaces and linear maps, we obtain the result. 

Programs written in the user language of kan/sml for algorithms in the 
present paper are available; please contact oaku@yokohama-cu.ac.jp. 



2 Computation of cohomology groups with 
coefficients in a locally constant sheaf of 
rank one 

A sheaf V on U is called a locally constant sheaf of rank m if for any x G U, 
there exists an open set W 3 x such that the restriction V\ w is a constant 
sheaf C^. 

Let fi, ■ ■ ■ , fd £ Q[ x ] be (not necessarily irreducible) factors of / satis- 
fying f — fi • • • fd- Let a±, . . . , aa be complex numbers which lie in a com- 
putable field. 

The left A n -module 

L(a) = Q[x,l/f]fr---f a d d 
is defined as follows; we define the action of and x^ by 

d(g/P 



d k .((g(x)/r)-m) = \^^j(g( x )/fP). m +^J-l. m , 
x k ■ ((g(x)/f p ) ■ m) = x k g(x)/f p ■ m 



where m = /j 21 • • • f^ d and g(x) is an arbitrary polynomial. In fact, we can 
easily check that 

d k ■ (x k (g/f p )m) = (d k x k ) ■ ((g/f p )m) 
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and hence our definition of the action is well-defined. 
The left A n -module 

P(a) = A,j? ■ ■ ■ f a / 

is the left A n -sub module of L(a) generated by m. 
Put 

V = Hom An {P{a),0?j n ), 

where O^J 1 is the sheaf of holomorphic functions on the complex manifold 
U = C n \ V(f). Here we endow U with the classical topology instead of the 
Zariski topology. When the left ^4 n -module P(a) is expressed as A n /I(a), we 
can regard V as a sheaf of holomorphic solutions on U of the system of linear 
partial differential equations 1(a); we have, for a simply connected open set 
u C lid, 

V{u) ~{fe 0^ n (u) \£f = for all £ E 1(a)} 
where the isomorphism is given by 

V(u) <p(l) G O a J(u). 

The C-vector space V(u) is one dimensional and spanned by the function 
fi 1 ' ' ' fd d since 1(a) contains 

d a .Mi. 

dk-Ys^T- (k = l,...,n). 

i=l fi 

Thus, V is the locally constant sheaf of rank one. 

Theorem 2.1 The cohomology group H k (U,V) is computable. 

This theorem is a generalization of Theorem [TT] . In fact, when a\ = 
. . . = ad = 0, the locally constant sheaf V is the constant sheaf Cjj. In 
order to prove this theorem, we need to generalize Procedure |1.3| to compute 
a left ideal 1(a) of A n such that L(a) = A n /I(a) where 1(a) is, intuitively 
speaking, the differential equations for (f^)fi 1 ■ ■ ■ f^ d with an appropriate 
nonnegative integer v. 

We introduce the Weyl algebra 

A-d+n = Q(£l, ■ ■ ■ ,td, d tl ,...,d td ,d 1 ,...,d n ). 
for our computation of 1(a). 
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Procedure 2.2 (Computing L(a).) 
Input: f,f 1 ,...,f d ,a 1 ,...,a d . 

Output: a left ideal 1(a) of A n such that L(a) = Q[x, l//]/" 1 • ■ ■ /d d ^ 
A n //(a). 

1. (Computation of the annihilating ideal of /f 1 • • • with indetermi- 
nates si, . . . , s d ) 

Compute 

(*j _ fj( x ) U = !> • • • > d )> + (i = 1, . . . ,n, j = 1, . . . , d)) 

n Q[ti<9 tl , • • • ,t d d td ](x,d x ). 

Replacing each tid ti by the indeterminate — Sj — 1 in generators of the 
intersection, we obtain the set 

G (-si - 1, ... , -Sd - 1) = ^a,, -s - 1), ... , Qh{x, d x , -s - 1)}. 



(Call Procedure |4.1| to compute the intersection of a left ideal and the 
subring Q[ti<9 tl , . . . , t d d td ] (x, d x ).) The left ideal I(s) of Q[ti<9 tl , . . . , t d d td ] (x, d x ) 
generated by Go (— s — 1) gives the annihilating ideal for /f 1 ■ ■ ■ f^ d . 

2. Compute 

(I(s),f 1 (x),...J d (x))nQ[s 1 ,...,s d \ 

by an elimination order x, d x > s±, . . . , s d . Let G±(s) be a set of gener- 
ators of the elimination ideal above. 

3. Choose a positive integer v such that the set 

(at - v, . . . , a d - v) - Z >0 (1, . . . , 1) 
is not contained in the zero set 

V(G 1 (s)) = {veC d \ g(v) = for all g(s) e G^s)}. 

4. Output 



1(a) :=G (-s-l) 



si^a±-v,...,s d ^a d -v 
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In the above procedure, 1(a) is the annihilating ideal of fi l U ■ ■ ■ fd d ~ u - 
The annihilating ideal of f? 1 ■ ■ ■ f% d can be computed as the ideal quotient 
1(a) : (A n f u ) through syzygy computation by means of Grobner basis. 

Let us present an algorithm to compute the co homology groups H k (U, V). 

Algorithm 2.3 (Computing the cohomology groups H k (U,V).) 



Input: f,fi,...,fd,ai,---,a d . 

Output: the cohomology groups H k (U,V). 



1. Call Procedure |27| with the input f,fi,---,fd, — a i, ■ ■ ■ , —Q-d- Get the 
output I (—a). 

2. Call Procedure |0|with the input / = /(—a). 



3 Computation of P(a) and its localization 

Put X = C n and let Y be an algebraic set of X defined by the polynomial 
/ G Q[x] with x = (xi, . . . , x n ). Let d = (d±, . . . , d n ) be the corresponding 
differentiations. We denote by Ox and V x = Ox(d±, . . . ,d n ) the sheaf of 
regular functions, and the sheaf of algebraic differential operators on X re- 
spectively (see, e.g., p. 15 and p. 70] and [0, p. 15]). We note that the set 
of the global sections T(X,T>x) coincides with C ®q A n , which is the Weyl 
algebra with coefficients in the complex numbers. We will denote it also by 
A n if there is no risk of confusion. 

In the sequel, we shall work in the category of algebraic Px- m odules and 
prove isomorphisms for sheaves of Px- m odules. Correctness of algorithms 
and procedures given in preceding sections follows by taking global section 
on X in isomorphisms of propositions. 

Put 

M=V(a)=V x fT---f a d d . 

The left coherent P^-module M. is a locally free 0x- m odule of rank one 
on X \ Y, which is called an integrable connection and -M\ X ^ Y has regular 
singularities along Y. As we have seen, in order to compute the cohomology 
group, we may compute j x -M [!//]• Our purpose in this section is to give 
a proof of correctness of Procedure |2.2j , which also gives an algorithm to 
compute the localization M. [1/f] '■— Ox[l/ f]<£>o x -M-. Ai[l/f] is a holonomic 
system on X (Theorem 1.3 of Kashiwara [|l(|) and coincides with M. on X\Y . 
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We outline a method to compute V(a)[l/ f] for given non-constant poly- 
nomials ft, ■ ■ ■ , fd an d a = (a\, ... ,0,4) with / := f\ - ■ ■ fd- Here, we assume 
that a,i lies in a computable field. 

Let s — (si, . . . , Sd) be commutative indeterminates and put 

C(s) := O x [8,l/f\f?~-f?, 

which we regard as a free Ox[s, l//]-module. Put V(s) := Px[s]/i si • • • fd d - 

Then the set of global sections T(X, C(s)) of C(s) coincides with C[x, s, 1/ f}}'] 1 ■ ■ ■ f d d , 

and that of V(s) with A^s}/" 1 ■■■f/. 



Definition 3.1 The (global) Bernstein-Sato ideal B(fi, . . . , fd) of Q[s] is 
defined by 



B(f h . . . , f d ) := {6(a) G Q[s] | 6( S )/f • • • ^ G • • • 



The step 2 of Procedure [2.2| gives an algorithm to compute the Bernstein- 
Sato ideal. 



Proposition 3.2 ( ||27|| ) There exist a finite number of linear forms Li(s), . . . , L K (s) 
in s with nonnegative integer coefficients, and nonzero univariate polynomials 
bx,...,b K , such that 

b(s) := bi(Lt(s)) ■ ■ ■ b K (L K (s)) G B(f x , . . . , f d ). 

In particular, for any a = (ai, . . . ,a d ) G C d , the intersection of {(ai — 
u, . . . , a d — v) I v G N} with 

V{B{f u f d )) := {s = (s u ...,s d )eC d \ b(s) = for any b G B(f u f d )} 
is a finite set. 

The following proposition tells us that if a is generic, then the localization 
C(a) of V(a) agrees with V(a). 

Proposition 3.3 Assume that a = {a±, . . . ,a d ) G C d satisfy that (ai — 
v,...,a,d — v) is not contained in V(B(fi, . . . , fd)) for any v = 1, 2, 3, . . .. 
Then V (a) = C(a) holds. In particular, the Ox-homomorphism f : V(a) — > 
V(a) is an isomorphism. 
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Proof: In the notation of Proposition |3.2| , there exist b(s) G B(fi, . . . , fa) 
and p(s) G A n [s] such that 

and b(ai — 1, . . . , a d — 1) ^ 0. Then we have 

A 01-1 • ■ ■ L ad_1 = 6(ai - !> • • " , 0.4 - lyM^i - 1, • ■ • , a d - ■ ■ • f a d *. 

Proceeding in the same way by using the assumption, we know that f\ x ~ v ■ ■ ■ fd d V 
is contained in V(a) for v — 1, 2, 3, . . .. This implies V(a) = £(a). [j 

Next, we shall see that the localization C(a) agrees with V(a\— Uq, . . . , a^— 
uq) for an integer vq determined by the zero set of the Berndstein-Sato ideal. 
In order to prove this fact, we need a lemma. 

Lemma 3.4 Ox[l/f) is a flat Ox-module. 

Proof: This should be well-known (e.g. this is a special case of Lemma 
1.1 of flTll). Here we give a direct proof. Let i : /C — > M be an arbitrary 
injective (9jf-homomorphism. Then for a section u of /C, we have 1 ® u = 
in K[l/f] = O x [l/f] ®o x K if and only if fu = for some v G N (cf. 
Lemma 7.2 of f24|). Hence we have 1 <g) u — if and only if 1 <g> i(u) = 0. 
This completes the proof. [\ 

Proposition 3.5 Fix an arbitrary a = (a±, . . . , a^) € C d . Le^ v$ be a positive 
integer such that {ax — v,...,a,i — v) is not contained in V(B(fi, . . . , /<j)) /or 
any integer v > vq. Then we have 

V{a)[l/f) = C{a) = V{a x -u ,...,a d - v ). 

Proof: Consider the short exact sequence 

— > V{a) C(a) — > C(a)/V(a) — ► 0, (3.1) 

where i is the inclusion. First note that (C(a) /V(a))[l / f] = 0. In fact, any 
section v of C(a) is written in the form v = gfi 1_u ■ ■ ■ fd d ~ u with g e Ox 
and v G N. Hence we have f u v G V(a). This implies (£(a)/7 : '(a))[l//] = 0. 
Since Ox\\/ f\ is a flat Cx-module, we have from ( ^.1|) an exact sequence 

— p(a)[l//]i^ £( fl )[l//] — ,0. 

Since £(a)[l//] = £(a), we have proved the first equality of the proposition. 
The second one follows from Proposition [O] since £(a) = C(a± — u , . . . , a<i — 
u ) (/ is invertible in C(a)). \\ 
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Proposition 3.6 Under the same assumption as Proposition |ff. j| , the T>x- 

homomorphism (specialization s = a) 

p : V(s)/(( Sl - at)V{s) + ... + (s d - a d )V{s)) — > V(a) 
is an isomorphism. 

Proof: Assume that a section u := p^fi 1 ■ ■ ■ f^ d of V(s) satisfies p(u) = 0, 
where u denotes the modulo class of u. Then there exist gx(s), . . . ,gd(s) G 
Ox[s\ and v G N such that 

«=E(^-^)/r-/r- 

3=1 

By the same argument as the proof of Proposition |3.3| , we can find b(s) G Q[s] 
and G such that 

Ks)fr u ---fr u = Q{s)K i ---r d d 

and 6(a) ^ 0. 

There exist ci(s), . . . , q(s) G C[s] which satisfy 

d 

6(a) — 6(s) = ^J(sj ~~ a j) c j( s )- 

3=1 

Hence we get 

= x>3 - pw» w/r" • • • + q,(a)pw/r • • • 

3=1 

= E( s 3 " (fc + A S1 • ■ ■ f/. 

3=1 

Since 6(a) 7^ by the assumption, we conclude that u G (si — a i )7 : '(s) + • • • + 
(s<2 — o>d)'P(s). Hence p is injective. The surjectivity is obvious. [J 

Let us consider the problem of finding the annihilating ideal of fl 1 ■ ■ ■ f^ d . 



20 



Let A d be the Weyl algebra on the variables t = (t 1 , . . . , t d ). We denote by 
A d V x '■= A d ® c V x the sheaf on X of the differential operators in variables 



(t, x) which are polynomials in t. We follow an argument of Malgrange (19 
for the case of d = 1. 

We can endow C(s) with a structure of left A^Px-niodule by 

t j (g^^)ft 1 ---f d d ) 
= g(x, Sl ,..., Sj + l,..., s d )f? ■ ■ ■ ft +1 • • • fT, (3.2) 

= -s j g(x,s 1 ,...,s j -l,...,s d )fl 1 ---fj"~---f d d (3.3) 
for g(x,s) e O x [s,l/f) and j = 1, . . . , d. 

Lemma 3.7 Let M be a sheaf of left ideals of A d T>x generated by 

tj-fi(x) {., 1 d). (3.4) 

d df- 

^ + E^ (i = l,-,n). (3.5) 

Then each stalk of M is a maximal left ideal. 
Proof: By a coordinate transformation 

t^tj-fii 30 ) C7 = 1 > •••><*)> 



x = x, 



we can reduce to the case where fi = . . . = f d = 0. In that case, the 
statement is obvious, [j 

Proposition 3.8 We have 

Af={ P eA d V x \pf^---f/ = 0}. 

Proof: It is easy to verify the inclusion C by using ( |3.2| ) and (^.3[). Since Af 
is maximal, we obtain the equality, [j 
We put 

X{s) := {p(s) G PxM I p(s)A ai ■ • ■ f d d = 0}. 
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Proposition 3.9 For a Zariski open set u of X , we have 



r(u,i(s)) 

= { P {- Sl - 1, . . . , - Sd - 1) I p(hd tl , . . . , t d d td ) e r(u, j\fnv x [t A , . . . , t d d td })}. 

Proof: By ( |3.4j) and (|3.5|) , we get the relations 

Sj = -d tj tj = -tjd t . - 1 {j = 1, . . . , d). 

Hence Px[s] is isomorphic to the subring V x [tid tl , . . . , t d d td ] of A d V X - This 
implies the conclusion. \\ 

Proposition 3.10 Procedure \2.2{ is correct. 

Proof: The correctness of step 1 follows from Proposition |3.9| . 

To verify the correctness of step 2 of Procedcure |2.2j , one has only to note 
that for b(s) G Q[s], we have b(s) G B(fi, . . . , f d ) if and only if b(s) — f 
belongs to F(X,I(s)). 

The correctness of steps 3 and 4 can be shown by taking global sections 
in sheaf isomorphisms given in Propositions |3.5| and p.6| . [] 

As to our experiments, it is more efficient that one eliminates d x first, 
and then elminates x in step 2 of Procedure 2.2 . However, even with this, 
the complexity of Procedure |2.2| is huge. 



4 Computation of the intersection of a left 
ideal and a subring 

In this section, we give a procedure to compute the intersection of the left 
ideal 

df- 

(tj - fj(x) (j = 1, . . . , d), -r^d t . + di (z = 1, . . . , n,j = 1, . . . , d)) 

in A d+n and the subring Q[tid tl , ■ ■ • , t d d td ](x, d x ) of A d+n . The intersection 
gives the annihilating ideal for f^ 1 ■ ■ ■ f^ d with the replacement tA l— > ~ s i ~ 
1. 

Procedure 4.1 Input: polynomials /i, . . . , in x — (xi, . . . , x n ). 
Output: a set of generators of the annihilating ideal X(s) of /* 1 ■ • • f d d . 
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1. Introducing indeterminates t = (t%, . . . , t^), u = (ui, . . . , Ud), v = 
(vi, . . . ,Vd), let I be the left ideal of A n+d [u,v] = Q[u,v](x,t,d x ,d t ) 
generated by 

tj-ujfj, (j = l,...,d), (4.1) 
d Bf 

d *i + Y, Q. u Ai 1 (4-2) 

1-UjVj, (j = l,...,rf). (4.3) 

2. Take any term order on A n+ d[u, v] for eliminating u, v. Let G be 
a Grobner basis of / with respect to this term order. Put Go = 
{P h ...,P k }:=GnA n+d . 

3. For each % — 1, . . . , k, there exist Qi G and vn, ...,!/«£ Z such 
that 

Si,un ' ' ' Sd t i/ id Pi = Qi{x, d x , tid tl , . . . , td,dt d ) 
holds, where Sj tU := d\. if u > 0, and Sj >v := tj 1 ' otherwise. Set 

G {s) := {Qi{x, d x ,s),..., Q k (x, d x , s)}. 
Output: Gq(— si — 1, ... , — Sd — 1) is a set of generators of X(s). 



Proposition 4.2 Procedure J^.l is correct. 



Proof: First, we must show that each element of Go can be written in the 
form as in the step 3 of Procedure O. Fix any j with I < j < d. Then 
the generators of / given in the step 1 are homogeneous with respect to the 
weight table Wj below: 



variables 


Xi,d Xi (l<i<n) 


tj 






V 3 


t k ,d tk ,u k ,v k {k^j) 


weight 





-1 


1 


-1 


1 






Moreover, the product of two operators preserves the homogeneity with re- 
spect to Wj. Hence each element of G is homogeneous with respect to Wj 
and free of u and v. This enables us to write Pj in the form as in the step 3. 

Now let us show that each Qi(x, d x , — s — 1) belongs to T{s) with the 
notation — s — 1 = (— Si — 1, ... , — s,i — 1). By the definition, Pj is contained 
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in the ideal generated by (|4.1|)-( fOD . Substituting 1 for every it, and v iy we 
know that Pj belongs to Af since it does not depend on u, v. Hence Qi(-s-l) 
belongs to X(s) in view of Proposition pT9| . 

Conversely, let p(— s — 1) be an arbitrary section of X(s). Multiplying by 
a polynomial, we may assume that p(tid tl , . . . , td<9t d ) belongs to the left ideal 



of A n+d generated by (|3.4j) and ( |3.5|) making use of Proposition |3]9| again. 
That is, there exist Rj, Si G A n+d such that 

d n I d o \ 

A , . . . , t^j = E ^ • - fi) + E s • K + E # (4-4) 

j=l i=l \ 3=1 aX « / 

We can homogenize the both sides of ( |4.4| ) by adding u with respect to the 
weight table Wj. By performing this procedure for every j = 1, . . . , d, we 
obtain a homogenization of (|4.4| ) with respect to all Wi, . . . , W^. The left 
hand side of this homogenization is in the form u 1 ^ 1 ■ ■ ■ u^ d p with nonnegative 
integers fix, . . . , [Aa since p itself is homogenous. Thus u^ 1 ■ ■ ■ u^ d p is contained 
in the ideal of A n [u] generated by ( |4.1| ) and ( |4.2j ). This implies that 

p — (1 — Ui 1 ■ ■ ■ U^V] 1 ■ ■ ■ v^ d )p + Ui 1 ■ ■ ■ U^Vi 1 ■ ■ ■ v% d p 

belongs to I. Since G is a Grobner basis of I with respect to a term order 
for eliminating u, v, there exist Ui, . . . , G A n+d such that 

k 

p{hd t ^ . . .,t d d td ) =J2 u i p i- 

i=l 

Since p and Pj are homogeneous with respect to each Wj, we may assume 
that so is Ui. Moreover, since the weight of p is zero with respect to each 
Wj, all Ui are written in the form 



Ui — Ui{tid tl} . . . } t d d td )S^ Vil ■ ■ ■ Sd,v id 

with some U[ G A n [tid tl , . . . ,t d d td ]. Hence p(s) belongs to the left ideal of 
A n [s] generated by G (s). This completes the proof. [| 

Example 4.1 Consider / = x Sl y S2 (l — x — y) S3 . X(s) is generated by 

ysi + ys 2 + ys 3 - s 2 - yxd x - y 2 d y + yd y , 
xsi + XS2 + XS3 — si — x 2 d x — yxdy + xd x , 
xs 2 + ys 2 + ys 3 - s 2 - yxd y - y 2 d y + yd y . 
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Note that this ideal is strictly larger than the ideal generated by trivial 
annihilators 

x(l-x-y)d x -x(l-x-y)(df/dx)/f, y(l-x-y)d y -y(l-x-y)(df/dy)/f. 



5 Twisted de Rham cohomology group 

In this section, we shall explain that computation of D-module theoretic 
integrals of L(a) gives the cohomology groups H k (U,V), which is nothing 
but what Grothendieck-Deligne comparison theorem says; the contents of this 
section should be well-known to specialists. However, they are not explicitly 
explained in literatures. 

First let us recall the integration functor for P-modules. In general, let A4 
be a left X>x- m odule (or, more generally, a complex of P^-modules) defined 
on X. Then integration of M. over X is defined by 



f M := RT(X,n x ®i x M) 



as an object of the derived category of C-vector spaces, where R and L de- 
notes the right and the left derived functors in the derived categories, V is 
the global section functor, and £l x is the sheaf of algebraic n- forms on X, 
which has a natural structure of the right Vx-modu\e and is isomorphic to 
Dx/(diDx + - ■ ■ + d n T>x) since X is the affine space. For i 6 Z, the z-th coho- 
mology of J x M. is denoted by J x -M, which is a C-vector space. R l T(X,J\f) 
is often denoted by H l (X,J\f). See, e.g., || and |T3j for an introduction to 
the mechanism of derived functors. 



Now put 



=1 fj 

Let M. be the left Px-module M. := Vx/1, where I is the left ideal generated 
by fdi — hi (i = 1, . . . ,n) with polynomials hi. Here, we note that hi satisfy 
the integrability condition 

( (hi/f) = Tr {h J /f) l<ij<n, (5.1) 



dxj dx 



and the function f® 1 ■ ■ ■ f% d is annihilated by the operators fdi — hi. M. has 
regular singularities along (the non-singular locus of) Y and also along the 
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hyperplane at infinity of the projective space P n . M. and V(a) are isomorphic 
as Px- m odules on X \ Y. In fact, both are simple holonomic systems and 
there exists a natural Dx-homomorphism of M. to V(a) which sends the 
modulo class of 1 G V x to f? 1 ■ ■ ■ f% d . However, these two modules are not 
isomorphic on X in general. 

The deRham complex DR(M[l/f\) of the localization M[l/f] := O x [l/ f}® 0x 
M. is defined by 

— Sl x ® 0x M[l/f] tt x ® 0x M[l/f] -±> • • • -±> ^® 0jt M[l//] — 0. 

(5.2) 

where d is given by 

n 

d(dx^ A • • • A rfx^ ® w) = fix,- A dx^ A • • • A rfx^ ® 

3=1 

for it G ,M[1//]. As Px[l//]-module (not as X>x- module!), there is an iso- 
morphism 

M[i/f] ~ ^[l/ji/^tv/K^ - ^r 1 ) + . . . + v x [i/f](d n - ^j- 1 )). 

Let P be a section of A4[l/f]. Then there exist Qj G X>x[l//] an d r G 
O x [l/f) such that 

n 

i=i 

Such r is determined uniquely. Then we define y^(P) = r. Hence 

<p:M[l/f\^O x [l/f\ 

defines an isomorphism as O x [l/ f]-modu\e. By transforming the complex 
( |5.2j ) by means of this </?, we get the following complex that is isomorphic to 

o — n° x [i/f] n^[i//] • • • nj[i//] — o, (5.3) 

where V, which is called the integrable connection, is defined by 



n ( du h \ 
Viudxi, A ■ ■ ■ A dxi, ) = — 1 — jr"U <ix 7 - A cfej, A ■ ■ ■ Adx 

pA dx 3 f J 

for it G C?x[l//]; in fact, we have 

_ „ du hj du 

d r u = udj + — = u-f + — 
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modulo PxtV/K^ - hf- 1 ) + . . . + V x [l/f](d n - hj- 1 ). Thus the integral 
fx f] = RT(X, ( |5.3|) ) is isomorphic to the complex 

o — > r(x ; n x [i/f}) r(x ; n^[i//]) • • • r(x ; nj[i//]) — o 

(5.4) 

since Q x [1 //] is a quasi-coherent Ox-module, X is affine and hence H k (X, Q p x [1//]) 



for k > (see, e.g., [|13| , p. 205, Proposition 1.2A, p. 215, Theorem 3.7] and 
|28f). The cohomology of this complex is nothing but the algebraic twisted 
de Rham cohomology with respect to the local system on X \ Y defined by 



the equation Vw = for u G Ox- When M = V(a) on X\Y, ( |5T4]) gives the 



algebraic twisted de Rham cohomology groups associated with the local sys- 
tem defined by V(—a), i.e. the cohomology groups of X\ Y with coefficients 
in the locally constant sheaf 

V := Hom Vx (V(-a), Of) = Hom Vx (O x ,Vf ® Vx V(a)), 

where Vf denotes the sheaf of holomorphic differential operators. In fact, 
by applying the functor Of®o x t° the complex (5.3), we obtain a complex 
of sheaves onX\7 whose A;-th cohomology group is 

{u e Of I Vu = 0} = Hom Vx (V(-a), Of) 

if k — and zero otherwise. 

The algebraic twisted de Rham cohomology coincides with the analytic 
one by virtue of the comparison theorem of Deligne @, p. 98 Theorem 6.2, 
p. 99 Corollary 6.3]. Let us summarize what we have explained. 

Theorem 5.1 (Comparison theorem, ||) 

H k (U, V) ~ H k ((T(X; n' x [l/f\), V)) ~ H k ~ n (X, DR(M[l/f])). 



As we will see in Proposition 3.1, we, moreover, have 



H k - n (X, DR(M[l/f})) ~ E k - n {A n l[d x A n + ■■■ + d n A n ) ® L An A n /{I{a))). 



Example 5.1 (Beta function) 

Putting X = C, we consider V(a) = T>xx ai (l — x) a2 for generic complex 
numbers a± and 02. We have V(a) = C(a) ~ T>x/ {p) with p = (x 2 — x)d x — 
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(ai + a^x + d\. The Bernstein-Sato ideal for x and 1 — x is generated by 
(si + l)(s2 + 1). The 6-function of the Fourier transform T>x/(p) with 

p = xd 2 x + (x + a x + a 2 + 2)d x + at + 1 



is s(s + ai + a2 + 1). Hence by applying Procedure [1.4| with fci = 0, we have 



— > + x^i) F /(F_! + x^i) 







and we get 



H°(U,V) =0, F 1 (?7,V) = C 

where {7 = {0, 1} and 

V(w) = {«g an H | du/dx = {-a x /x + a 2 /(l - x))u] 
for a simply connected open set w. Note that 



H\U,V) 



C 



x 



i 



' x(l-x) 



vc 



~ C 



X 1 — X 



dx 



where V = d+ (ai/x — 02/(1 — x))dxA. The beta function should be regarded 



as 



x ai (l -x) a2 <p 



where ^ = ^fy e H\U,V). 



Example 5.2 For generic complex numbers a%, . . . , a m , we consider V(a) = 
T>xIViLi(x — Ci) ai where Ci,...,c m are distinct points in C. By applying 
our algorithm, we can see that H l (U,V) = C m_1 and H°(U,V) = where 
U = C\ {ci, . . . , c m } and V = Hom Vu (V(-a), Off 1 ). See § for details on 
these cohomology groups and hypergeometric functions. 



Example 5.3 (Counting the number of bounded chambers by P-module 
algorithms) 

We consider a collection of hyperplanes 



Li(x) — ^ c ij x j + c i0 
3=1 



0, (i = l,...,m) 
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in R n and put / = HiLiLi(x). For complex numbers ai, . . . ,a m , we con- 
sider V(a) = T> x Y\™ =1 L i (x) a \ The number of bounded chambers in U = 
R n \ UZi{x | Li(x) = 0} is equal to the Euler number of H k (U,V) (see @, 
p. 47 Theorem 2.13.1] and @). Although there are several algorithms in com- 
putational geometry to count the number, this number can also be counted 
by our purely algebraic algorithm. Table [I] is an example of computation of 
Euler numbers by our algorithm and implementation. 



Table 1: dim c ^(X \ Y, C) 



/ 


i = 2 


% = 1 


i = 


Euler ch. 


xy 


1 


2 


1 





xy(x + y + 1) 


3 


3 


1 


1 


xy 










■(x + y + 1) 


6 


4 


1 


3 


■(x-y-2) 











where X = C 2 = {(x,y)}, Y — {f — 0}. 

Of course, our method is far from efficient. However, it is rather surpris- 
ing that purely algebraic computations in the ring of differential operators 
can evaluate the number of bounded chambers in a given hyperplane arrange- 
ment. 



6 Computation of integration 

Let Ai be holonomic Dx-niodule defined on X := C n . In this section, we 
explain a method to translate the computation of integrals p x Ai to that of 
the restriction H l ((T> x /(xiV x + - • ■+x n T> x ) ®A n M), where Ai is the Fourier 
transform of Ai. Our discussion together with the algorithm of computing 
the restriction in p5] prove the correctness of Procedure [O] and consequently 



the correctness of steps 2, 3 and 4 of Algorithm |1.2| . 

Let us denote by VL X the sheaf of regular (algebraic) i-forms on X. We 
use the notation d = (di, . . . , d n ) with di := d/dxi. Let us denote by DR(Ai) 
the complex 

— > Q° x ® 0x M V x ® Qx M ^ Q x ® Qx Ai — ► 0, 
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where d is defined by 

n 

d(dXk 1 A ■ • • A dx ki ® u) = ^ dxj A dx kl A ■ ■ ■ A dx ki ® {dju) 

3=1 

for u G A4. Here we regard Q l ®o x -M as being placed at degree i — n. In 
particular, the cohomology groups of DR{V X ) are given by 



W*(DR(Z> X )) 



fix if z = 0, 
otherwise. 



Hence we have an isomorphism 

^lx®v x M ~ DR(V x )®v x M = DR(M). 

Since Q x ®o x -M is a quasi- coherent Ojf-module and X is afflne, we have 

H k {X^ x ® 0x M) = (fc>0). 

Hence by using the standard argument for the sheaf cohomology, the integral 
is explicitly represented by a complex j x A4 = RT(X; DR(.M)), which is 
equivalent to 

— ► (A Z n )® z M (A Z")® Z M -^-> U (A Z n )® z M — ► 0, (6.1) 

where M := .M) and 

d(e i:1 A . . . A e ik ®u) — e 7 - A e it A . . . A e ih ® (djw) 

i=i 

with the unit vectors ei, . . . , e n of Z n . 

The Weyl algebra A n has a ring automorphism $ defined by 

= ${di)=Xi (i = l,...,n), 

This $ naturally defines a new left A n -module M := $(M), which is called 
the Fourier transform of M. Since M. is holonomic, M belongs to the Bern- 
stein class of A n -modules (cf. [[5], p. 125]). Since the Bernstein class is invariant 
under the Fourier transform, we know that M. := T> x ®A n $(M) is a holo- 
nomic Px-module on X. By applying $ to the complex (5.1), we obtain 
another complex 

— ► (A Z n )® z $(M) (A Z n )O z $(M) U (A Z n )<g> z $(M) — ► 0, 

(6.2) 
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A e 



(XjU) 



where 

n 

5(e il A . . . A e ik ® u) — ^2 ej A e ix A 

i=i 

Since the complexes (|6.1| ) and ( |6.2|) are isomorphic, we have only to compute 
the cohomology groups of ( |6.2| ). Here note that fl6.2| ) is a complex defining the 
restriction of Ai to the origin of X. Thus, we have the following proposition. 

Proposition 6.1 We have for any i, 

H\X, DR(M)) ~ W((A n /( Xl A n + ...+ x n A n )) ® L An M). 

Note that Ai is specializable to the origin (i.e., a nonzero 6-function 
exists) since Ai is holonomic (cf. |]18|). Hence the cohomology groups of 
) are computable by steps 2 - 6 of Procedure L4 as shown in [25]. Thus 



each P x Ai is computable as a finite dimensional vector space and we obtain 
the following proposition. 



Proposition 6.2 Procedure 1.4 is correct. 



The heart of Procedure [L4| is the truncation of the resolution with re- 
spect to a filtration defined by the weight vector w by a root of 6-function 
[25]. Let us briefly explain the idea. Let Ai be a holonomic T>- module and 
b(s; x) be the 6-function (or indicial polynomial) along x\ = 0. Then, xy : 
grk+i(Ai) p —> grk(M) p is bijective if b(k;p) 7^ 0. Here gr(Ai) is the graded 
module associated with the weight vector w = (—1, 0, . . . , 0; 1, 0, . . . , 0). 
Hence, in order to obtain the kernel and the image of the map Xi, we may 
truncate the high degree part and the low degree part of the filtration of M. 
with respect to the weight vector w. In order to obtain all the cohomology 
groups of the restriction, we need a diagram chase to determine the degree 

2|, section 5] and |25| . 



of the truncation. As to details, see 



By Proposition |3.10| , Proposition \L.2\ Theorem |5.1| and Proposition |6.2| , 
we obtain the following theorem and complete our proof of Theorems and 



Theorem 6.3 Algorithms l.i and 2.S are correct. 



We close this section with the following theorem, which generalizes The- 



orem |2.1| under the condition that the coefficient sheaf is expressed in terms 
of A n -module M. Note that we did not require this condition in Theorem 
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Theorem 6.4 Let M be an A n -module (A n ) p /I where I is a left submodule 
of {A n ) p . We assume that A4 = T>xM is regular holonomic on X and also 
regular along P n \C n and that the singular locus of M. on X is given by f — 
with a polynomial f . The cohomology group H k (U, HoniD i7 (C[/, T>fj i ®A n M)) 
is computable where U = C n \ V(f) . 

Proof: An algorithm to compute M[l/f] is given in Tensor product and 
localization] (see also [pi] , Propositions 7.1 and 7.5]). The rest of the proof 
is same as that of Theorems |0| and |2j] because M.\. , can be regarded as a 
regular connection ([]ISL Theorem 2.3.2]). \\ 

Example 6.1 If M. is not regular holonomic on P n , then the comparison 
theorem no longer holds. For example, put 

M = V x /(d x + 2x), X = U = C. 

The operator d x + 2x is not regular at x = oo. We can see that 

H°{X,DR(M)) = C, H-\X,DR(M)) = 

by applying our integration algorithm. Now, take y? G Homx> x (Ox, -M an ). 
We may assume that / = <^(1) belongs to (9 an since d x = — 2x in Ai 31 * 1 . We 
have d x f = in Ai an , which means that fd x + /' G (d x + 2x). Then, we have 
/'// = 2x and hence / = v?(l) = ce x G Ai an for a constant c. Therefore, we 
have Hom Vx (Ox, -M an ) ~ C. On the other hand, 

H\X, C) = + H°(X, DR(M)) and H°(X, C) = C ^ H'^X, DR(M)). 

7 Computation of cohomology groups on the 
complement of an algebraic set when its 
algebraic local cohomology group vanishes 
except for one degree 

The purpose of this section is to establish a connection between the de Rham 
cohomology of C n with an algebraic set removed, and the integration of 
modules over the Weyl algebra. We use the algebraic local cohomology groups 
lying in between these two objects. The contents of this section except the 
last theorem should be well-known to specialists. 
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Put X be an n-dimensional non-singular algebraic variety over C and 
let Y be an arbitrary algebraic set of X. The algebraic local cohomology 
group H[y](Ox) (in the sense of Grothendieck) is the i-ih derived functor 
of the functor F^y] of taking the support (in the algebraic sense). This is a 
holonomic T>x-module (Theorem 1.4 of ||16|| ). When X = C n , algorithms for 
computing the algebraic local cohomology groups have been given in for 
the case where Y is of codimension one and f34|, [2|| for the general case. 

Proposition 7.1 Let Cx be the constant sheaf on X with stalk C. Then 
there is an isomorphism 



RT Y (X,Cx) ^ [ RT [Y] (O x )[-n] 



where [— n] denotes the shift operator in the derived category. In particular, 
if TilyAOx) = for i d, then for any i e Z, there is an isomorphism 

ri—n — d 

H Y {X-C X )^ / nf Y] {o x ). 

Proof: The algebraic and the analytic de Rham complexes are defined by 

DR(Ox) := n x ®^ x Ox, 
DR an (C^ n ) := ^<g)£ M £> an 



x i 



where fi^f denotes the sheaf of holomorphic n-forms. Then we have 

J x RT [Y] {O x ) = RT{X-n x ® L Vx RT [Y] {Ox)) 

= RT(X;RT [Y] (DR(Ox))) 
= RT Y (X;DR(Ox)). 

On the other hand, there are two distinguished triangles and a morphism 
between them: 

RT Y (X;DR(O x )) — > RT(X;DR(O x )) RT(X \Y;DR(O x )) 

Rr Y (X;DR an (0^)) — ► RF(X;DR an (Of)) -±i> RF(X \ Y; DR an (OJf )). 

Here the vertical homomorphisms except the leftmost one are isomorphisms 
by virtue of the comparison theorem of Grothendieck [12]. Hence the leftmost 
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vertical homomorphism is also an isomorphism. Moreover the complex de 
Rham lemma implies DR an ((9^ 1 ) = Cx[n]. Consequently, we get 

RT Y {X;DR(O x )) = RT Y (X;DR™(O a «)) = RT Y {X; C x )[n\. 

This completes the proof. \\ 

From the above proposition and the isomorphism H l (X\Y; C) ~ H^^X; C) 
(see, e.g., [H| p. 212, exercises 2.3]), we obtain 

Corollary 7.2 Assume 7i l ^(O x ) = for i ^ d. Then for any i > 1, we 

have an isomorphism 

i—n—d+l 



H\X\Y-C X )~ f Hf Y] {O x ). 



If Y is non-singular, we can also relate the de Rham cohomology of X \ Y 
to that of Y itself: 

Corollary 7.3 Assume that Y is non-singular and of codimension d. Then, 
for any i > 1, there exists an isomorphism 

H l {X \ Y; C x ) ~ H l+1 - 2d {Y- C Y ). 

Hence, H l (Y; Cy) is computable for any i > 0. 

Proof: Let i : Y — > X be the embedding. Then by the Kashiwara equivalence 
(cf. Ul4] , p. 34, Theorem 1.6.1] and we have an isomorphism 7~tf Y ] {O x ) = 
l + O y . Thus by using Proposition [TT| , we obtain 

n—d 



rt—n—a 

h y (x-c x ) ~ / nf Y] {o x 



i—n—d 

l+O y 

X 

—n—d 

o Y 

Y 

~ W-^iY-Cv). 

Combining this with the preceding corollary, we are done. [| 

In |[34| , U. Walther gave an algorithm to compute the local cohomol- 
ogy groups 7i Y {M) with a Cech complex under the condition that M. is 
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(fi • • • /d)-saturated. Since Ox satisfies this condition, we can compute al- 
gebraic local cohomology groups TClyAOx) for any i > where Y := {fi = 
. . . = fd = 0} . Another approach to compute algebraic local cohomology 
groups of M. with a resolution and without the condition of saturation is 
given in f25fl . Thus, we have two algorithms for the next theorem. 



Theorem 7.4 The cohomology groups H l {X \ Y; Cx) for any i > is com- 
putable if Y = V(fi, . . . , fa), fi € Q[xi, . . . , x n ] and if HfyAOx) vanishes 
except for one j . 

This theorem generalizes Theorem [TTTJ under the condition on vanishing 
of the local cohomology groups Ti^iOx)- Note that if d = 1, then this 
condition always holds. 
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